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ABSTRACT 


The relationship between the bending moment and the 
through-plane shear force in the vicinity of a mechanical fas- 
tener at failure was determined. Experiments were conducted 
on 4-inch wide flat plate aluminum and graphite-epoxy composite 
speciments that modeled portions of a wing skin along a spar 
and along a rib. The composite specimens were either 8-ply 
or l6-ply balanced layups and were simply supported at two 
opposing edges and free along the other two edges. The fas- 
teners were pulled normal to the plates, and the maximum force 
at failure was measured for specimen lengths varying from two 
to six inches between supports. 

The aluminum plates failed by formation of a plastic hinge 
across their width and showed little sensitivity to through- 
plane shear. The 8-ply spar specimens cracked across their 
width and also were relatively insensitive to through-plane 
shear. However, failures of the rib specimens were confined 
to a region near the fastener, where the fastener pulled 
through the plate, and showed much greater sensitivity to 
through-plane shear. 

Two analyses were made; one for small elastic deflections 
of a thin orthotropic plate, and another for a beam in the 
elastic range. A mesh generator for a finite element model of 
the plate around the fastener was also developed for the 


computer program ADINA. 
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І. INTRODUCTION 


А. BACKGROUND 

In recent years considerable effort has been devoted to 
the development of new materials for the purpose of improving 
aircraft performance. One family of these materials is fiber 
reinforced composites. Many studies have shown that use of 
advanced composite materials in aircraft structures can result 
in significant weight savings due to their low density, high 
modulus character. Furthermore, the anisotropic nature of 
fiber reinforced composites challenges the aircraft designer 
to exploit their directional nature in order to realize even 
greater savings. 

Graphite-epoxy composite materials have recently been speci- 
fied for major structural elements in two U.S. Navy and Marine 
aircraft currently in advanced deveiopment stages.  Specifi- 
cally, the wing skin, leading edge extension, trailing edge 
flap, rudder, horizontal stabilizer, and vertical fin are to 
be constructed of graphite-epoxy in the F-18, [1]. The under- 
lying wing substructure is to be constructed of metal in a 
traditional spar-rib arrangement. The second aircraft, the 
AV-8B, is the advanced version of the AV-8 Harrier. This 
aircraft will have the entire wing box, consisting of the skin 
and the torque box substructure, constructed of graphite- 


epoxy, [2] 
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Mechanical fasteners are specified to attach the wing skins 
to the underlying wing substructure for both aircraft to allow 
removal of the skins,[3]. 

Extensive tests under various conditions have been con- 
ducted to assess the risks associated with graphite-epoxy 
materials. For example, strength degradation due to moisture 
absorption, high temperature environments, galvanic corrosion, 
and fatigue have been determined for the service environment 
of both aircraft, [3]. These tests have shown that design 
strain levels for both the Harrier and the F-18 wing designs 
provide adequate safety margins for the detrimental effects 
of elevated temperature and moisture absorption to be encoun- 
tered in service. 

ШОЧ ОТОП ЕО ene ability of the composite structure to 
withstand design loads under a wide range of service conditions, 
some attention has been given to the vulnerability of air- 
craft composite structures in a combat environment. Tests 
have been conducted to study the response of advanced composites 
to ballistic impact damage and the concomitant reduction in 
strength. A discussion of these tests is given in Reference [4]. 

Much data are available and a great deal is known about 
the behavior of metals at failure. Theories of failure have 
been tested time and again and failure models have proven va- 
luable in predicting stress levels that lead to failure. With 
the advent of composite materials, attempts to model the behav- 


ior of composites and the failure mechanisms have led to some 
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understanding of failure modes. However, significant gaps 
exist in the literature. Of particular interest here is the 
failure of a composite structure at a metal fastener in a 


particuiar mode called fastener pull-through. 


B. FASTENER PULL-THROUGH 

Consider a composite plate en fastened to the under- 
lying structure with mechanical fasteners, as in the F-18 and 
Harrier wing designs. Primary considerations in selecting skin 
thicknesses include in-plane shear, tensile strength and criti- 
cal buckling stresses. Joints are primarily designed to trans- 
fer in-plane tensile and compressive loads in the skin to the 
underlying spar-stiffener-rib structure. 


Three different failure modes at a joint are shown in Figure 


l. These are bearing, shear-out, and tensile failure. These 
failure modes are discussed in detail in Reference [5]. A 
fourth type of joint failure is illustrated in Figure 2. Іп 


this failure mode the fastener is pulled through the plate by 
a large through-plane shear force. Although very little atten- 
tion has been given to this type of failure in metals, this 
failure mechanism may be significant in fibrous composites 
under certain loading conditions. 

For example, through-plane shear stress exists in the 
Skin of pressurized fuel tanks. Strain levels are usually 
small under normal design conditions, and the through-plane 


shear stresses at the joint are insignificant. However, if the 
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aircraft takes a hit in combat, circumstances may arise which 

subject joints around the fuel tanks to unusually high through- 
plane shear stresses. The extent to which the joint is able to 
carry these stresses may have significant impact on the vulner- 


ability of the aircraft. 


S HYDRAULIC RAM 

Consider the problem of an integral wing fuel tank which is 
partially full of fluid when impacted by a high speed projectile. 
The kinetic energy of the projectile is transferred to the fluid 
in the tank as it passes through, [6]. The result of this en- 
ergy transfer has been characterized by three phases. The 
first phase is the propagation of a shock through the fluid. 
The second is the formation of a vapor filled cavity behind the 
projectile. The third phase is an oscillatory phase during 
which the cavity expands and contracts as the system returns 
to an equilibrium state. These phenomena, illustrated schema- 
tically in Figure 3, have been referred to as the hydraulic 
ram effect. 

The intense fluid pressure created by the hydraulic ram 
has a two-fold effect on the tank wall. First, the pressure, 
acting perpendicular to the wall, creates a moment distribu- 
tion through the plate. Secondly, the through-plane shear 
forces that arise must be reacted by the fasteners. This 


Situation is shown in Figure 4. 
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Reference [7] describes hydraulic ram tests that were 
conducted on fluid containing aluminum tanks and one manifesta- 
tion of the failure of the tank walls was an "unzipping" of 
fasteners. Although the fasteners that were used in these 
tests were not typical aircraft fasteners, questions were 
raised about such connection failures. If aluminum plates 
and rivets failed in this through-plane mode, then it may be 
possible that composite plates will fail in the same way. 
Although composites are generally pound-for-pound stronger than 
aluminum in-plane, they may be considerably weaker in the 
through-plane direction. In addition to the bending moment 
and through-plane shear induced by the hydraulic ram loading, 
aircraft skins will be in tension, compression, and in-plane 
Shear due to maneuvering when the hit occurs. This combina- 


tion of stresses can lead to a significantly lower failure load. 


ЕФ OBJECTIVES 

Several aluminum and graphite-epoxy plates have been 
tested to determine the sensitivity of failure to the com- 
bination of through-plane shear and bending. moment due to 
hydraulic ram loading. Specifically, the objectives were to 
determine the force required to pull a fastener through the 
plates for a range of values of the bending moment at the 
fastener. A failure curve relating pull force and bending 
moment similar to a yield curve for limit analysis was sought. 


Since reinforcing fibers do not run in the through-plane 
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direction in composite plates, it was expected that through- 
plane shear strength would be low compared to aluminum plates. 
Macroscopic failure mechanisms in both materials were 
observed. In order to verify experimental observations, two 
analyses were made. One analysis was for small, elastic de- 
flections to determine the magnitude of the errors in moment 
distribution associated with treating the specimens as beams 
instead of plates. Another model applies to beams in the 
plastic range. Additionally, a mesh generator for use in 


the elastic-plastic finite element code ADINA was developed. 
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ТТ. EXPERIMENTAL APPROACH AND SET-UP 


A. THE MODEL 

From Figure 4 it is apparent that the hydraulic ram 
loading on a plate creates a complex internal stress dis- 

e ution throughout the plate and particularly at the 
fastener. To further complicate matters, hydraulic ram is a 
dynamic problem involving large pressures applied over a 
short time interval. In order to facilitate the experimental 
and analytical approach, some simplifying assumptions were 
made. First, the dynamic effects of the hydraulic ram were 
not addressed, although it is known that thev are of great 
importance in inducing catastrophic failures at fuel tank 
boundaries. Secondly, it was assumed that there would Бе по 
in-plane stresses present. Again, it was recognized that these 
will play an important role in joint failures. These assump- 
tions are justified here on the basis that this is a first 
study in an ongoing program. Future studies will consider 
these effects. 

Figures 5A and 5B illustrate the modeling of the skin 
with the test plates. The underlying structure may be envi- 
sioned as either a rib, Figure 5A, or a spar, Figure 5B. The 
difference has importance only with respect to the composite 
plates. For example, the zero degree ply direction was as- 


sumed to run parallel to the main wing spars, [1]. 





Figure 5C shows the loading set-up to simulate the com- 
bination of bending moment and shearing force that exist at 


the fastener due to the hydraulic ram loading. The edges at 


y t b/2 were assumed to be free and those at x = O and 

х = a were assumed to be simply supported. Тһе free condition 
was assumed because, near the ribs and spars, moments due to 
the hydraulic ram loading would be smaller in the y direction 


compared to the moments in the X direction. 


ЕЕ TEST EQUIPMENT 

In order to load a specimen as suggested by Figure 5C, 
a devide was constructed in the Naval Postgraduate School 
Machine Shop. The simple support condition was provided by 
triangular knife edges fashioned from steel bar stock. A 
means was provided by which the supports coulá be positioned 
from One to five inches from the center of the specimen in 
one-inch increments. By moving the supports outward, the 
moment at the fastener becomes larger for a A pull torce: 
Thus, supports close together provide a relatively larger 
pull force and smaller moment whereas supports farther apart 
give a relatively smaller pull force and larger moment. The 
supports were attached to a rigid steel base to facilitate 
attachment of the device to the moveable head of a Riehle 
200,000 lb. testing machine. The device is shown in Figure 6 
with an aluminum plate in place for a test pull. The fasteners 


were passed through a hole drilled in the specimens and at- 
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tached to a Clevis device which in turn was attached to a 
Baldwin SR-4 5000 lb. capacity load cell. 
The load cell was connected to a locally constructed 


wheatstone bridge. A power supply and digital voltmeter 


were connected to the bridge circuit in the usual way and cali- 


brated to read an integer multiple of the load applied at 

the fastener. Bridge output was also monitored on a strip- 
chart recorder calibrated to read pounds of pull. The strip- 
chart recorder was most useful in providing a time history of 
the applied loading and assisted in recording significant 
events which occurred during the loading. Some aluminum 
specimens were fitted with strain gauges, which were wired to 
a wheatstone bridge switching device. Gauge output was moni- 
tored by a second digital voltmeter calibrated to read micro- 
inches per inch. An additional channel of the strip-chart 
recorder was available to monitor the strain gauge output 
selected at the manual switching station. Figure 7 shows the 


entire experimental facility. 


ER TEST SPECIMENS 
l. Aluminum Specimens 
The aluminum specimens were machined from sheet stock 
to a four-inch nominal width.  Lengths of the specimens were 
varied from seven to fourteen inches, depending on the in- 
tended location of the supports. A hole was machined in the 


center of the specimen to provide a nominal two-thousandths 
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inch interference fit for the fastener as specified in 
Reference [8]. 

The fasteners that were used for the tests are the 
Hi-Tigue Hiloks that have been specified as skin fasteners 
for portions of the F-18. These HLT-328-8-18 fasteners 
are often specified for applications which require high fatigue 
resistance. Their unique Hi-Tigue feature, however, was not 
germane to the tests. Figure 8 shows the dimensions of the 
fastener. 

2. Composite Specimens 

The composite plate specimens were fabricated by the 
author from Hercules AS-3501-6 twelve-inch pre-preg tape. 

The tape was a nominal ten mil thickness with low resin con- 
tent. Five l6-inch x l6-inch plates were constructed using 
полен о-ршу [0/45/90]. or 16-ply [0/+45/0,/+45/90] „ balanced 
layups and were prepared in accordance with Reference [9]. 

The composite preparation equipment is shown in Figure 9. 

Four-inch wide test specimens were cut from the 16- 
inch square plates after post-curing. Some 8-ply specimens 
were prepared such that the 0-degree ply was perpendicular to 
their axis (the spar specimens). The remaining 8-ply and the 
l6-ply specimens had the O-degree ply parallel to the axis 
(the rib specimens). The plates were drilled and fitted with 
the same fasteners as the aluminum specimens, although inter- 


ference fitting is not specified for this type of material. 
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Room temperature material properties from Reference [5] 
are tabulated along with computed flexural rigidities in 


Table I. 
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III. ANALYTICAL AND COMPUTER MODELS 


A. ANALYTICAL PLATE MODEL 

The test specimens are idealized as thin plates with 
two edges simply supported and two edges free, as shown in 
Figure 10. The deflections and moment distribution can be 
approximated by assuming the load to be uniformly distributed 
over the shaded rectangle of area uV. The governing equation 
for a balanced, symmetric laminated plate is given by [6]. 


D w 1) 


+ + 
a ОС 2 


where W is the plate deflection, 4 is the applied pressure 


load and subscripts denote partial derivatives. D , D 


, D , 
11 ша 292 


and D , are stiffness coefficients. In the special case of 
isotropy we have, 
= р = р; р з ур; р  - D(l-u)/2 (2) 
11 242 12 66 
Тһе solution to Equation (1) for the portion of the plate 


prst shown in Figure 10 is given by [9]. 


w= > t (a, PA cosh (mmy /a) + (3) 


MID © o 
+8, (Mr y/a)sinh (Mr y/a)} sin (mr x/a) 
where, 


a, = 4Pa4(-1) (п-1) /2 sa 1) (4) 


5 5 
T m D Uv 


and a is the span of the plate. 
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Considering the unloaded portion of the plate beyond 


line ts, assume a deflection surface of the form, 


w! «УҒА! cosh(my/a) « B', (m /a) sinh(mmy/a) (5) 


3 c! Sinh» тула) + D'.mmy/a )cosh(nmy/a)) sinmmx/A) 


The six constants A_, B., AP DE, must satisfy the 


m 
free edge boundary conditions at y = b/2, and continuity 
Eng tine tS. Applying continuity conditions at y = b/2 


gives 


i 1 1 
W= We.: W = W . W =h - W =W 
f ч и“ 222 к" ‘ххх ‘ххх (6) 


It is shown in Reference [9] that 
А ЗА! _ 
m m = a (Ym sinh 2Y mn cosh 2Y mm) (7) 


= 1 == 
pu B m 27122 (cosh 2 т) 





- 2 
= a, (y„ cosh 2Y n sinh 2Y mm) 
E | 
Do = a sinh 2ү,)/2 
where 
ym = MTY (8) 
4а 


Two more equations are available from the boundary conditions 
at y=b/2 in order to solve for the six constants. Along the 
edge y-b/2 the moment м, and the shearing force E must 


vanish. Hence 


- -(D w, Б у=. о 59) 
22 [2 


М, | 
Y y=b/2 yy 


= - (D 


W, +(D +2 D м, = 0 
кс КО D. н н.) 


2 6 6 XXV 
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b titution of Equation (4) into Equation (9) gives, 

А! колы: 

= ——— 
m пт") чук 

(1-9 Ja, cosha,-(o,-5)sinha, ED sinha, -2 cosha, 


(o, -1)sinha, (1-9 )cosha,, 


е 1 = . E 1 А . 
С m o )sinha, D E (1 0 Ja cosha,+2 sinha, ) 








x 
ЕГІ = E й 1 Е m 
8 х! o )cosha, D m (1 o JayQSinha, (о, 3) cosha,, ) 
й 
where 
к = (1-0 ) (1-0 294 * [2 (p, -1) + (o -1) (0, -5)] cosha, sinha, 
. mmb Ly = i жээ? (Dr2+2D66). Е Юг 
m шт | Баро 8 D22 ds Dii 


The quantity that is of particular interest here is the 
Ex er bution of M along the line x з а/2. A plot of M. /M 
versus 2y/b is shown in Figure 11 for the aluminum specimen, 
and each of the three types of composite specimens for various 


values of a, where 
tu 
4 (11) 


М = 


PLASTIC LIMIT ANALYSIS FOR A BEAM 
If the aluminum specimen is considered as a beam, then 


the total bending moment is given by Reference [12] as 


M “ил 20(х) 42 012) 
-t/2 
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Where b is the width and t is the thickness of the beam. 
Making the usual assumptions that plane sections remain 
plane, the strain at any point is given by 

€ = KZ (13) 
Where Kk is the curvature of the middle surface, z= 0. For 
moments below a certain critical value, Mu all stresses 
are elastic and 

= Ee =EKZ (14) 


Substituting into Equation (12) gives 
Бек (15) 


If the moment is increased until at the outer fibers, 
2 = +t/2, is equal to the yield и ‚ the maximum 
elastic moment is given by 

Me = poto, (16) 
Since only the outermost fibers of the beam are yielded, the 
beam will continue to carry additional load through the cen- 
tral fibers until they too are yielded. The situation is 
shown schematically in Figure 12. The limiting moment is 


Case (D) of Figure 12 where the entire cross section is 


yielded and, 


M = с. Бъ: (17) 


The moment in the beam is largest at the center of the 


s anderes value 1S ра/4. Substituting Equation (11) into 
Equation (17), the maximum pull that the beam can carry, Ру, 
is given by р ас 2 ¿2 (18) 

y уа 
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Coe Gon GENERATOR FOR A FINITE ELEMENT MODEL 

The theory of matrix structural analysis is discussed at 
length in Reference [13], and a discussion of various three- 
dimensional elements can be found in Reference [14]. The 
studies presented in Reference [14] by Clough demonstrate that 
hexahedral elements with nodes on the sides adequately repre- 
sent the bending of a simply supported plate. Clough used 
a Single layer of 20 node elements and achieved the exact 
deflections. He showed that more complicated elements, such 
as curved tetrahedra, could perform as well, but that the 
formulation time was significantly greater. Further, it was 
shown that, while two-point Gauss quadrature rules gave exact 
results for rectangular prisms, three- or four-point rules 
were required for skewed elements in order to achieve the 
desired accuracy. 

In this study, 16 node isoparametric elements were chosen 
to model specimens. The objective of the finite element 
effort was primarily to develop a mesh generator which could 
later be used to analyze the laminated anisotropic plates. 

The finite element code used was ADINA, developed by 
Bathe [15] in 1975. ADINA incorporates a library of four 
basic elements and allows selection of any combinationof 15 
material models. The library material model that seemed 
best suited was the isotropic elastic-plastic model with strain 


hardening.  Tensile tests were conducted on aluminum coupons 
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of three different alloys in order to accurately provide 
proper input data for ADINA. The results of the tensile 
tests were used throughout the investigation as baseline 
material properties, and the stress strain curves for three 
different aluminum alloys are shown in Figures 13A through 
13C. The bilinear approximation to the stress strain be- 
havior of the materials is indicated by the broken lines. 
Symmetry of the test specimens allowed some simplifi- 
cation of the model. The analysis includes one-quarter of the 
plate. The finite element model of the specimen is shown in 
Figure 14. The loads can be applied as concentrated loads at 
nodes along lines (A) and (B) of Figure 14. At x = 0, re- 
straints can be applied such that no displacement in the x 
direction is permitted. А+ у = 0, по displacement in the y 
direction is permitted. At x = a, the simple support condition 
dictates that no displacement in the z direction is permitted. 
The only other constraints involve the nodes at z= 0 and on 
the first ring of elements closest to the fastener. The assump- 
tion is made that the fastener is essentially rigid with respect 
to the aluminum. Consequently, no displacement in either the 
X or y directions is permitted for these nodes. 
A description of the mesh generator, and its use is in- 
cluded in Appendix A. Present limitations internal to the 
ADINA source code preclude the possibility of obtaining a 


successful nonlinear analysis due to array dimension. 
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IV. RESULTS AND DISCUSSION 


A.  ALUMINUM SPECIMENS 

All of the aluminum specimens tested failed by forming a 
plastic hinge across the center of the span, as shown in 
Figure 15. None of the fasteners pulled through, and there was 
no visual indication that the through-plane shearing stresses 
had any effect on the failure mode. 

Table II contains the ultimate moment and pull force 
based upon the plastic limit analysis and the actual maximum 
moments, Me, and pull forces, Pg, obtained during the experi- 
ments. Note that due to the fact that a - 4 inches, Pe and 
Me have the same numerical value. The ultimate moments ob- 
tained experimentally are seen to be somewhat higher than 
those predicted by the plastic limit analysis of section III. 
Part of the difference may be due to the fact that the theory 
applies to elastic-perfectly plastic materials, whereas the 
materials tested (particularly the AL7075T6) exhibited some 
degree of strain hardening. Further, the limiting moments 
predicted by the theory are average moments, whereas the 
moment is not distributed uniformly across the beam, as shown 
in Figure ll. If the ultimate pull force is taken to be that 
force required to raise the moment M along the Line X = 4/2 
to the ultimate moment, then the predicted values of pull, 
Pe, more closely resemble those obtained experimentally. 


In any case, the tests indicate that these aluminum plates fail 
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at the joint in a plastic hinge mode. This hinge effect has 
previously been observed in hydraulic ram experiments, (71. 

Figures 16 and 17 show a plot of the stress versus applied 
load for two of the locations which were provided with strain 
gauges on the 6061 aluminum specimen. The stress is normal- 
ized by the yield stress for the material, and the load is 
normalized by the theoretical Pe for the material and support 
conditions. The figures also contain the theoretical stress 
for a beam with a concentrated load of the same magnitude, 
and for the flat plate. For these tests, the simple supports 
were 4 inches apart. It can be seen that at the point x = 2, 
y = 1.040, the beam approximation does not predict stress levels 
accurately. The plate theory is considerably better. Further 
away from the centerline of the plate, at X = 2.735, Y =.735, 
the plate theory overestimates the experimentally obtained 
stress levels. This suggests the possibility that formulation 
of the plate problem may not be accurate near the fastener. 
The assumption that the fastener head behaves like a uniform 
load acting over a rectangular area may not be refined enough 
to accurately predict the stresses in the vicinity of the 


fastener. 


B. COMPOSITE SPECIMENS 
The results of all tests on the composite specimens are 


presented in Table III. M was determined for each of the 


ült 


five plates by applying a line load across the center of a 
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four-inch wide by seven-inch long specimen cut from each plate. 
The line load was applied by a steel cylinder attached to 

the Riehle testing machine, as shown in Figure 184. The fas- 
Bener passes through a hole in the cylinder and the hole in 
m plate and attaches to the Clevis. Thus, pulling the fas- 
tener causes the cylinder to bear against the plate along the 
EN uot the cylinder. This eliminates the concentration of 
emr ough-plane shear at the fastener. Р 1+ was determined 

for each 16-inch square plate by pulling a fastener through a 
two-inch by two-inch specimen restrained by a two-inch by six- 
inch by five-eighths-inch thick steel plate with a circular 
hole eighty-thousandths of an inch larger than the radius of 


the fastener head. The А? set-up is shown in Figure 18B. 


Je 

Two complete sets of experiments were run on the 8-ply rib 
specimens in order to establish that the results were not at- 
tributable to defects in the quality of the plates from which 
the specimens were cut. The mean difference was less than 7 
percent. 

БЕЗЕУ Ра plot of the experimental data in nondimen- 
sional form. The abscissa is the moment, M, normalized оу 
Mult The ordinate is the pull force at failure, "en normal- 
ied by Р 1167 The extent ог mode of failure is indicated in 
the figure. 

Јаше са Gonsist of photographs of the 8-ply spar, 
8-ply rib, and 16-ply rib failed specimens. The sequence of 


photographs for each group of specimens is from small moment 





mM mostly through-plane shear, to small through-plane shear 
and mostly moment loading conditions. The extent of the failure 
zone is indicative of the sensitivity of the specimen to through- 
plane shear. Note that the failure zone for specimens loaded 
with relatively small moment is primarily confined to a re- 
gion around the fastener hole. As the moment was increased 
(by lengthening the span) this failure zone extended in the 
КОШ стоп until it covered the entire width and resembled the 
failure of the iine loaded specimen. 

Figure Z1A is a plot of the data from a different point of 
view. The abscissa in this plot is the theoretical elastic 
local moment at the fastener head, 02 normalized by Malt 
Росте 218 is a plot of the theoretical elastic moment at the 


edge of the plate at the center d the span, normalized 


Hedge’ 
by Mult! The theoretical values were obtained from Figures 


E cand 11B at Y = 0.4 inches. In both cases the ordinate is 


the pull force at Failure, Pe, normalized by Р 1167 
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V. CONCLUSIONS 


A.  ALUMINUM SPECIMENS 

Ша пас пов росвір1е to pull the HL-328-8 fasteners through 
the aluminum test plates. The failure of the aluminum joints 
was manifested entirely by the formation of a plastic hinge 
and the resultant inability of the specimen to carry further 
load. 

There is essentially no comparison between the way the 
aluminum specimens and the composite specimens failed when 
 естей to the combined effects of moment and through-plane 
shear. Failure of the aluminum plates was insensitive to the 
through-plane shear force induced by the head of the fastener. 

The assumption that the plates behaved like beams when sub- 
jected to a force on the fastener does not appear to be valid, 
with the exception of the special case where moment (or support 
Spacing) is very large. It was shown from plate theory that 
the effect of the concentrated load upon the moment distribution 
is more significant where support spacings are small. The ap- 
proach taken here, therefore, requires some refinement before 
it can accurately predict stress levels. However, use of a 
finite element model should aliow designers to more accurately 
predict stress levels that lead to failure. A wide variety 
of loading conditions and fastener geometries may be studied 


En this manner. 








COMPOSITE SPECIMENS 

It is evident from Figure 20 that the weakly-oriented 
Bemposite spar specimens, 52 and 53, were largely insensitive 
to through-plane shear. The location of the failure points on 
Пише 10 is also indicative of a lack of sensitivity since 
these specimens failed at a load that was very near the ulti- 
mate moment. However, the rib specimen showed a significant 
sensitivity to through-plane shear. Thethickest specimens 
showed the greatest degree of sensitivity, and this result may 
be cause for concern. 

The manner in which pull-through failure occurs is complex. 
It is suspected that the ultimate moment is attained over an 
area near the fastener, causing failure of the specimen in 
that vicinity. However, the distribution of moments over the 
centerline is such that a portion of the plate outside this 
failure zone maintains sufficient strength to carry further 
load. The through-plane shear forces must still be carried by 
the composite around the fastener. Thus, as the load is in- 
creased further, the fastener pulls through the specimen. In 
those specimens which were loaded with large moment and small 
ШИГ through, the distribution of moments is flatter along the 
Ме тле, запсе the supports are further apart. That is, the 
effect of the concentrated load induced by the fastener head 
is less pronounced. Conseuqnetly, once local failure around 


the fastener occurs, only a small increase in pull force is 
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Pa 


I med tO increase the moment to its critical value along 
the entire width of the specimen. Thicker and wider speci- 
mens, in which the ultimate moment is larger, may show greater 
sensitivity to through-plane forces, and therefore be less 
likely to fail across their entire width. It would be in 

(и Stiff materials where "unzipping," if it is to occur, 
would be most likely. Current wing design concepts for por- 
mens of the F-18 wing call for thicknesses from .3 to about 

Wm пслез. Clearly, as thickness gets larger, the through- 
plane shear strength gets larger and there is a point where the 
actual through-plane shear loads can be carried by the com- 


pesite. 
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VI. RECOMMENDATIONS 


A. ANALYTICAL METHODS 

A mesh generator has been developed to model the simply 
supported- simply supported-free-free plate loaded over a 
ШИЛ” region Under the fastener head. This has been done 
using an elementary loading condition and an isotropic material. 
A great deal of understanding could be gained if the program 
were modified to include orthotropic elements. The current 
program allows generation of up to 16 layers of 16-node iso- 
tropic isoparametric elements. The node generation and ele- 
ment connectivity would remain unchanged. 

A further sophistication could be realized by ‘treating the 
nodes near the fastener as fixed and applying pressure loads 
as predicted by the computer program Satans [16. Satans has 
the ability to predict pressure distributions on a plate that 
arise in a fluid as a consequence of hydraulic ram from a 
Seecttic threat projectile. Other modifications could include: 
1. The ability to change the geometry of the fastener head. 
пе ability to take into account the flexibility of the 


fastener head. 


Са 


ono In plane tensile, Compression, and shear loads 


Characteristic of maneuvering flight. 


ES E PERIMENTAL PROCEDURES 
No data were obtained at the lower moment ranges because 


шим го 7 CG move the support arms of the machine closer 
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together than two inches. MI Edcceom derectoshould be cor- 


rected by affixing the support arms to a base which permits a 


mI Trance of motion, and does not restrict selection of dis- 


tances between arms to one inch increments. 


Investigation of the sensitivity of composites to a more 


complicated stress state which includes in-plane stresses is 


required. Further study is also required on factors which 


affect sensitivity, such as: 


lg 


Fastener head geometry variations, such as counter-sunk 
heads. 

Fastener-hole buffer materials and seals, such as O-rings. 
Degree of interference fit or tolerance between fastener 
and hole. 

ptrecte of'high cycle fatigue. 

Existing delaminations or other damage. 


Effects of softener strips. 
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APPENDIX A 


A. GENERAL 

The mesh generator is extremely simple to use. Where possi- 
ble, variable names are exactly the same as in Reference [15]. 
In addition to the source deck, only four data cards are re- 
quired to generate a data deck that can be used as input for 
ADINA. By simply removing the load cards, the deck can also be 
used with the preprocessor program PSAP1l. The generator is 
specific and therefore limited in that it currently has the 
ability to generate nodes along 5 radial lines in one quadrant 
of a plate that has a hole in the center. Other nodes are less 
restricted. The group of elements closest to the fastener hole 
are fixed in size by fastener parameters that are input. Тһе 
size of the rest of the elements is determined by the number of 
divisions desired along the line that lies at the center of the 
span of the plate. The capability currently exists to generate 
up to 16 layers of 16 node isotropic isoparametric elements. 
However, the source program ADINA requires minor internal modi- 
fication to accommodate the large volume of storage required by 
a non-linear analysis. A program is currently underway at the 
Naval Postgraduate School to provide ADINA with this increased 


сарара 1 лЕу. 


Br REOUIRED INPUT FOR MESH GENERATOR 


Input data for the generator requires only four cards. 
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Input Card Format (1615) 


COLUMN 

Ша, Дед 
oreet 
Il S15 
oO 2 0 
21485525 
GD 
SUS З 
Do ea 
dun d ye 
46 - 50 
SS 55 
56 = 60 
ові UG 
DONO 
DE D 
6 S0 


VARIABLE DESCRIPTION 
трови) 
DIESE) 
IDOF (3) 
IDOF (4) 
TBO (Ss) 
IDOF (6) 
NEGNL 
MODEX 
NSTE 
TPRI 
NPB 


LAYERS Number of layers of elements 
desired. 


NEGL 
NLOAD 
NLCUR 


NPTM 


Input Card 2 FORMAT (20A4) 
1 - 80 


БИН ALPHANUMERIC TITLE 


DESIRED ON 
ADINA AND PSAPI PRINTOUTS. 
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INPUT CARD 3 FORMAT (8F10.0) 


l - 10 PROP (1) Young's Modulus 

11 - 20 PROP (2) Poisson's Ratio 

21 - 30 PROP (3) Simple tension 

31 - 40 РКОР (4) Strain Hardening Modulus 
A Length of 

4l - 50 half span between simple 
supports 

5l - 60 ВІ Radius of Fastener 

6l - 70 PMAX Maximum load desired 

71 - 80 PINC Desired load increment for 

load data cards. Different 


from DT which is load step 
internal to ADINA, and deals 
with the method of solution. 


INPUT CARD 4 FORMAT (8F10.0) 


иг! 0 R Басе сі Hastencr Head 
11 - 20 T Thickness of specimen 
21 - 30 DT Loadstep increment 
Allow to default (0.0) 
3l - 40 TSTART Allow to default (0.0) 
41 - 50 FAC Allow to default (0.0) 
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<k KMESH GENERATOR PROGRAM LISTINGH****%% 
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